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We investigate the QCD ghost model of dark energy in the framework of Brans-Dicke cosmology. 
First, we study the non-interacting ghost dark energy in a flat Brans-Dicke theory. In this case we 
obtain the EoS and the deceleration parameters and a differential equation governing the evolution 
of ghost energy density. Interestingly enough, we find that the EoS parameter of the non-interacting 
ghost dark energy can cross the phantom line (wd = — 1) provided the parameters of the model are 
chosen suitably. Then, we generalize the study to the interacting ghost dark energy in both flat and 
non-flat Brans-Dicke framework and find out that the transition of wd to phantom regime can be 
more easily achieved for than when resort to the Einstein field equations is made. 



I. INTRODUCTION 



o : 

<N ■ 

Many observational evidences such as data from type la supernovae [l[ , cosmic microwave background (CMB) 
and SDSS 0, Q have convinced people to accept that our universe is currently experiencing a phase of accelerated 
Q\ , expansion. Based on the Einstein's theory of gravity we should introduce a new type of energy with a negative pressure 
to provide such an epoch of evolution. This unknown energy component which push the universe to accelerate usually 
' (— i ' called dark energy (DE) in the literature. Of course, the simplest candidate for DE is the cosmological constant 
originally introduced by Einstein to explain the static behavior of the universe at that time. However, cosmological 
constant suffers coincidence and fine tuning problems. Also a dynamic approach to dark energy has been proposed 
, avoiding such problems and in favor of experimental evidences seeking a time vary ing equation of state (EoSVparameter 
u>d = Pd/pd- Examples of such models a re q uintessence @, 0], K-essence (8l-[Toj. phantom [ll|, quintom [33], tachyon 
holographic [13] and agegraphic DE |l5j . 

■ In recent years there has been a new attention to the so called "scalar-tensor gravity". Scalar-tensor models of 
| gravity are those introduce a scalar field modifying the Einstein's theory of gravity. These models are mainly seen to 

retrieve at the low energy limit of string theory. However, the history of these models back to many years ago. The 
scalar tensor theories have started with the work by P. Jordan in 1950 16]. A prototype of such models was proposed 
! by Brans and Dicke in 1961 [13]. Their aim for presenting this model was to modify Einstein's theory in a way 
admitting the so called "Mach's principle". To this end they proposed a new scalar degree of freedom to incorporate 
the Mach's principle into general relativity [I?} • The Brans-Dicke (BD) theory can pass the experimental tests from 

■ the solar system jl8j and provide an explanation for the acceleration of the universe expansion [l9| . 

Most DE models introduce new degrees of freedom in addition to those exist in the standard model of cosmology. 
Introducing such degrees of freedom need investigation about their nature and new consequences in the universe. 
Hence, it seems so impressive and economic if one can explain the DE puzzle using already presented fluids and fields 
in the literature. Ghost dark energy (GDE) is an example of these models which uses the so called Veneziano ghost 
field to explain the recent acceleration of the universe [2(| Hl|. Ori gina lly the so called "Veneziano ghost field" was 
presented as a solution to U(l) problem in effective low energy QCD |22| - [25| . Taking into account the ghost field leads 
some consequences in the vacuum energy density in a dynamic spacetimc or a spacetime with nontrivial topology, 
while such a field seems to be un-physical in the Minkowski spacetime and has no contribution to the vacuum energy 
density. In curved spacetime the ghost field gives rise to a small vacuum energy density proportional to Aq CD H, 
where H is the Hubble parameter and Aq CD is QCD mass scale [2l|. It is shown that such a vacuum energy density 
is capable to drive a phase of acceleration and can be considered as a dynamical cosmological constant [26|, |27j ■ It 
is important to note that Freidmann equation for GDE model in Einstein gravity is similar to the self-accelerating 
branch of the Dvali-Gabadadze-Porrati (DGP) braneworld 28, 29], that is both of them have the form H 2 oc aH+pw, 
[2ll |30| . However we would like to emphasize here that in DGP braneworld the term proportinal to H appears due 
to the inclusion of large extra dimension in the gravity theory and hence modification of the Einstein field equations 
which yields an extra degree of freedom in the theory, while GDE model is totally embedded in standard model and 
general relativity, one needs not to introduce any new parameter, new degree of freedom or to modify gravity. As 
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a result the origin of the term aH in GDE model completely differs from previous models such as self-accelerating 
DGP braneworld scenario [30j . 

Although, it is a general belief that the current curvature of the universe is negligible and mostly the universe is 
considered with a flat geometry, recent observations support the possibility of a non-flat universe and detect a small 
deviation from k = [3l| . For example evidences from CMB and also supernova measurements of the cubic correction 
to the luminosity distance favor a positively curved universe [U HH . In addition, some exact analysis of the WMAP 
data reveals the possibility of a closed universe [Hj|. 

Various aspects of GDE have recently investigated. A generalization of GDE model in a non-flat universe was 
discussed in (35J . Tachyon and quintessence reconstruction of GDE model were studied in [H, H3| • Since GDE model 
belongs to a dynamical cosmological constant, it is more natural to study it in the framework of BD theory than in 
Einstein gravity. In this paper, we study a cosmological model of late acceleration based on the GDE model in the 
framework of non-flat BD cosmology. 

This paper is organized as follows. In the next section, we present the GDE model in the flat BD theory. Interacting 
GDE model in a flat BD theory is discussed in section III. In section HVj we generalize the study to the universe with 
spatial curvature. We summarize our results in section Ivl 



II. GHOST DARK ENERGY IN BD THEORY 

The action of BD theory, in the canonical form, can be written [38| 

S = J d 4 x^g (- ±-q> 2 R + l -g^d^d v 4> + L M \ (1) 

where R is the scalar curvature and <f> is the BD scalar field. The non-minimal coupling term <f> 2 R replaces with the 
Einstein-Hilbert term R/G in such a way that G~$ = 2-k(/) 2 /uj, where G e ff is the effective gravitational constant as 
long as the dynamical scalar field 4> varies slowly. In the Jordan frame, the matter minimally couples to the metric 
and there is no interaction between the scalar field <f> and the matter fields. The equations of motion for the metric 
g^v and the BD scalar field are 

4>G\w = -8nT^ - ^ \ 4>#<l>,v - ^g^,x^ x ) - <t>;w + (2) 
2wT3" 



2- 

u< t> = 7^;T^\ (3) 



where TM, is the matter energy-momentum tensor. Our aim in this paper is to consider the GDE in the Friedmann- 
Robertson- Walker (FRW) universe which is described by the line element 



ds 2 = dt 2 _ a 2 [t) + r 2 dn 2 ^j , (4) 

where a(t) is the scale factor, and k is the curvature parameter with k = —1,0, 1 corresponding to open, flat, and 
closed universes, respectively. A closed universe with a small positive curvature (f2& ~ 0.01) is compatible with 
observations [3{|- Using metric (fj|, the field equations ^ and ((31) reduce to 

i> 2 (4 + " 2 4) -i"** - - 5 (' + =) * - (6) 

* + 3//*-i-f° + H 2 + i)* = 0, (7) 
2uj \a a z J 

where H = a/a is the Hubble parameter, po and pd are, respectively, the energy density and pressure of DE, and 
Pm is the pressureless dark matter (DM) density. 

Consider the FRW universe filled with DE and pressureless matter which evolves according to their conservation 
laws 

p D + 3Hp D (l + w D ) =0, (8) 
Pm + 3Hp M = 0, (9) 
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where wjj is the EoS parameter of DE. In this section we want to consider the GDE in a spatially flat spacetime in 
the BD framework. The ghost energy density is proportional to the Hubble parameter [2lJ 



Pb = aH. 



(10) 



Here a is a constant of order Aq CD where Aqcd ~ lOOMeV is QCD mass scale. With H ~ 10 33 eV, Aq CD £T gives 
the right order of magnitude ~ (3 x 10~ 3 eU) 4 for the observed DE density [2l| . 

To determine the evolution of the universe filled by the pressureless matter and GDE, through the equations ([5][7]), 
we still have another degree of freedom in analyzing the set of equations. Based on the previous experiences in the 
BD theory it is known that usually the BD scalar field cj) has a power law relation as 



4> = (f> a(t) £ 



(11) 



A case of particular interest is that when e is small whereas ui is high so that the product su> results of order unity 
[4Cj |. This is interesting because local astronomical experiments set a very high lower bound on ui (4lj |: in particular, 
the Cassini experiment implies that u > 10 4 [H, Taking the derivative with respect to time of relation (fTT|) we 
obtain 



^ = £- = eH. 

6 a 



Using Eqs. (JTTJ) and (Ti"2"]) . the first Friedmann equation §5§ becomes 



H 2 (l - — e 2 + 2e) + — = —^{p D + p M ). 
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We introduce the fractional energy densities corresponding to each energy component as usual 
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where we have defined 



3fH 2 



Using (|10p we can rewrite Eq. p6[) as 



tin = 



AUJOL 

3<j) 2 H' 



Based on these definitions, Eq. p3|) can be rewritten as 

7 = n D + n M - n k , 

where 



7 = 1 



2lo 



(12) 

(13) 

(14) 
(15) 
(16) 

(17) 

(18) 

(19) 
(20) 



Taking the time derivative of relation (1101) , and using the Friedmann equation (|13[) as well as the continuity equations 
© and ©, we find 



tH = * = -H 
PD H 



3 3£Idwd 
£+ 2 + 2 — 



Substituting this relation in continuity equation © we obtain the EoS parameter of GDE, namely 



w D 



7 



2j-n D 



2 ±-l 

3 



(21) 



(22) 
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(gj = 10000,Q d = . 72;Q k = 0) 



0.01 0.02 0.03 0.04 0.05 



. 5 
0.25 



■0.25 
-0.5 
■0.75 





















0.2 0.4 0.6 


0.8 1 







FIG. 1: The first figure indicates wd against e in the flat non-interacting GDE. In the second figure the deceleration parameter 
(solid) and w e g (dashed) are plotted versus Qd for e = 0.002. 



It is worthy to note that in the limiting case e = (cj — > oo) we have 7 = 1 and hence the BD scalar field becomes 
trivial; as a result Eq. (j2"2"]l reduces to its respective expression in flat standard cosmology 27] 
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w D = 



2-0, 



(23) 



The solar-system experiments give the result for the value of u is w > 40000 [18] . However, when probing the larger 
scales, the limit obtained will be weaker than this result. In Ref. [42j], the authors found that oj is smaller than 40000 
on the cosmological scales. Also, Wu and Chen [43| obtained the observational constraints on BD model in a flat 
universe with cosmological constant and cold DM using the latest WMAP and SDSS data. They found that within 
2cr range, the value of u> satisfies u> < —120.0 or u > 97.8 j4^|. They also obtained the constraint on the rate of change 
of G at present 



G 

- 1.75 x KT^yr -1 < — < 1.05 x lO -12 ^ -1 

G 

at 2cr confidence level. So in our case with assumption (fTT|) we get 

G 



G 



= - = eH < 1.05 x 10 



■12 , 



yr 



(24) 



(25) 



This relation can be used to put an upper bound on e. Assuming the present value of the Hubble parameter to be 
Hq ~ 0.7 we obtain 



e < 0.01. 



(26) 



The GDE model in BD framework has an interesting feature compared to the GDE model in Einstein's gravity. 
It was shown that in standard cosmology based on Einstein's theory, the EoS parameter of the noninteracting GDE 
cannot cross the phantom line Wjj = — 1 and at the late time where Vlp — > 1 approaches — l[27l|. However, in the 
BD framework, requiring condition wd < —1 leads to 7(^ + 1) < £!d- Choosing ilo = 0.72 for the present time, 
this inequality valid provided we take e = 0.002 which is consistent with observations. This indicates that one can 
generate a phantom-like EoS for the noninteracting GDE in the BD framework. 

In addition to the EoS parameter of the GDE we can also study the effective EoS parameter, w e s, which is defined 

as 



Pt Pd 

W c tf = = • 

Pt Pd + Pm 



(27) 



where pt and Pt are, respectively, the total energy density and total pressure of the universe. As usual we assumed 
the dark matter is in the form pressureless fluid (Pm=0). Using relation (fT9|) for the flat case one can find 



n, 



-w D 



2 7 - Oi 



Hf-i 
3 



(28) 
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It is also interesting to study the behavior of the deceleration parameter defined as 



1-4- (29) 



Inserting from (f2~Tj) into (|29|) yields 



1 3 D io D 



Substituting Eq. (|2"21 in the above relation one can easily reach 

i m D / •>: 



2 2(2 7 - fi/j) V 3 



- 1 . (31) 



Let us study some special cases of interest for the deceleration parameter q. If we take fip = 0.72 for the present 
time and choosing e = 0.002 and u> = 10 4 we obtain q = —0.36, which is consistent with the present value of the 
deceleration parameter obtained in [44|. This choice of parameters lead Wr> = —0.79 and w e s — —0.58. The evolution 
of the deceleration parameter q and w c g are plotted in the second part of Fig([T]). A close look to this figure reveals 
that the universe enters acceleration phase when Cljj = 0.49. The effective EoS parameter at flu = 0.49 becomes 
Weft — —0.33 while wd = —0.66. We can also study the behavior of the deceleration parameter in the early stage of 
the universe where -^C 1 . In this epoch q = h + e which indicates that the universe was experiencing a phase of 
deceleration at the early stage of its evolution due to the domination of the DM component. 

Once again we can see that in the limiting case e = (7 = 1) the above relation restores the deceleration parameter 
of GDE in Einstein's gravity [13] • 

1 3 n D 

q= 2 "2(2^)- (32) 

Finally, we obtain a differential equation governing the evolution of GDE from early deceleration to late time accel- 
eration. To do this we take time derivative of the relation (TTgj) and use Eq. (f2T)|) . The result is 

Ct D =CL D H(l + q-2e). (33) 

Inserting q from (l3Tj) and using relation il' D = H , we obtain 

i-n D 



-2s 



(34) 



where prime denotes the derivation with respect to x — In a. In the limiting case e = (7 = 1) one recovers the result 
obtained in 12711. 



III. INTERACTING GHOST DARK ENERGY IN FLAT BD THEORY 



Most of the models study the dark side of the universe consider the evolution of DE and DM separately. This means 
that pdm and pde are separately conserved. However, recently there has been a lot of interest in interacting models 
of DE, since observations detects a signal of interaction between DE and DM. For instance, observational evidences 
provided by the galaxy cluster Abell A586 supports the interaction between DE and DM [45( ■ Beside there is no any 
reason against the interacting behavior of DE and DM and many authors has explained Lagrangians leaded to an 
interacting approach [46| . Any conservation in the physics can be explained by a symmetry in the Lagrangian level 
and there does not exist such a symmetry through the known Lagragians explaining the dark side of the universe. 
These points convince us to discuss a version of GDE in which there exist an interaction between DE and DM. To 
this end we can write the conservation equation for different components as 

PD + 3H PD (l + w D ) = -Q, (35) 
pM + 3Hp M = Q, (36) 

where Q denotes the interaction term and we take it as Q — 3b H(pm + Pd) with b 2 is a coupling constant. Such a 
choice for interacting term implies the the DE and DM component do not conserve separately while the total density 
is still conserved through 



p + 3H{p + P) = 0, 



(37) 
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(ai= 10000, Q D = 0. 72; Q k = 0) (e = 0.002,b=0.15) 
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FIG. 2: In the first figure wd is plotted against e. Solid line corresponds to b — 0.05 and dashed one for b = 0.2. In the second 
figure the deceleration parameter (solid) and the w e g (dashed) are plotted versus JId- These figures are plotted in the flat, 
interacting GDE case. 



where p = po + Pm and P = Pd- In this section we would like to discuss the interacting GDE in a flat universe. To 
this goal, once again we take time derivative from (|10p to get (|21l) . Next, we obtain the L.H.S. of Eq. (1211) from the 
conservation equation and the R.H.S. from 
the EoS parameter 



T3| and also taking into account Q = 3b 2 H(pM + Pd) and (| 1Q[) we obtain 



w D 



7 



2rf-Sl D 



2e 

y 



2jb 2 



(38) 



In the limiting case e = (u — > oo) or 7 = 1, the BD theory reduces to the Einstein gravity and hence Eq. 
recovers its respective expression in standard cosmology (27j 



1 



w D = 

The effective EoS parameter can be written as 



2-n 



D 



"(I 



9., 



7 



-w D 



ftD~' 



I"' 



2 7 b 2 

n 



D 



In this case, using the method of the previous section, the deceleration parameter is obtained as 



D 



2(2 1 -n D ) 



2s 



- 1 



276 1 



(39) 



(40) 



(41) 



In order to obtain an insight about various features of the interacting case, we try a suitable choice of parameters. 
Taking Op = 0.72, b = 0.15 and e = 0.002, we obtain q — —0.41 for the present value of the deceleration parameter 
which is in good agreement with recent observational results fjij . Also this choice of the parameters result in 
wd — —0.84 and w e s — —0.62 for present time of the universe and one can see the consistency of the latter with current 
observations. Transition from deceleration to acceleration take places at firj = 0.46 and at this time u> e g = —0.33 
while Wd — —0.71. One may take a look on Fig. [5] for a better insight into the model. The first part of this figure 
indicates that one can generate a phantom-like behavior provided e < 0.01. Further, we can have a close look at the 
behavior of wd at the late time where Qd ~^ 1- In this limit and using the same values of e and b we find wo = —1.06, 
indicating that the interacting GDE in the BD framework can cross the phantom line in the future. In this limit with 
a same set of parameters w e s = —1.09 which leads a super acceleration and may result in a big rip as fate of the 
universe. 

The equation of motion of GDE is obtained as 



D 



3 -(l « 



2jtf 



- 2e 



(42) 
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(a)=10000,Q D = 0.72;Q k = 0.01) (e = 0.002,b=0.35) 
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FIG. 3: In this figure wd is plotted against e. Solid line corresponds to b = 0.05 and dashed one for b = 0.2. The second 
figure shows the evolution of deceleration parameter (solid) and also w c r (dashed) versus Qd- These figures are plotted in the 
presence of the curvature as well as the interaction between DE and DM for GDE in BD framework. 



IV. INTERACTING GHOST DARK ENERGY IN NON-FLAT UNIVERSE 

Historically, one of the problems leaded to pioneering theory of inflation was the so-called "flatness problem", 
flatness problem was explained by assuming a very huge rate of expansion in the first stages of the cosmic evolution 
namely "inflation" . After the inflation it is generally believed that the universe is spatially flat, that is the inflation 
practically washes out the effect of curvature in the early stages of cosmic evolution. However, in the past decade 
several observational evidences have been observed in contrast to the flatness assumption of the universe. In the 
context of inflation it is argued that the flatness is not a necessary consequence of inflation if the number of e-foldings 
is not very large [47j . Besides, the parameter fl k represents the contribution to the total energy density from the 
spatial curvature and it is constrained as —0.0175 < fife < 0.0085 with 95% confidence level by current observations 
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48]. In addition, there are several observational evidences which supports a closed geometry for our universe in 
-|34|. Based on these implications it seems interesting to investigate the different DE models while the curvature 



term is present. So, in this section our main aim is to discuss a non-flat FRW background in the presence of the 
interacting GDE in BD theory. Taking the curvature into account, the Friedmann equation is rewritten as 

JH 2 + A = i^L ( fim + pD ) . (43) 
a z 3<p z 

Following the method of the previous section we obtain different parameters in the presence of the curvature term. 
Taking the time derivative of the Friedmann equation (|43[) and using (|19|) . we find 



H 2 _ 7 7 



3 3 Qdwd 
£ + - + -■ 



2 2 7 +ftfc. 

Combining the above equation with Eq. (|3"5]) . after using Eq. (|2"Tj) . we obtain the EoS parameter as 



(44) 



»»--5-V +->+£<•'+«'>)• (45) 

27 - il D \ 37 3 7 Q D J 

One can easily check that the above equation reduces to the EoS parameter of non-flat interacting GDE in the 
Einstein's framework provided e — (7 = 1). 

Using the same method as the previous sections, w e g in a non-flat background reads 

w cff = — — ^r w D- (46) 
7 + "fc 

This relation can be obtained using Eq. (j2"7| as well as (|T9"f in the non-flat case. Replacing (14"51) in the above relation, 
we get 

w eff = -J7y n w . a , 1 - IT - -ri 1 + ) + 77- (7 + ■ (47 

(27 - il D )h + "fc) V 37 3 7 S1 D / 
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The deceleration parameter can be calculated by substituting Eqs. (|44|) and (|45|) into (|2Q|) . We find 



(1 



7 ' 



27£ 



1 



2 n D -27 



2(2? -£)d) 



3 



Also the equation of motion of GDE can be obtained by replacing Eq. 



D 



G?ln a 



1 + (1 + ^) 

7 



2 7 e 



27 



3fi 



D 



2(2j-Q 1 



in (|33J}. The result is 

ft* 2b 2 
3t" 



+ rj-(7 + ^fe) 



2c 



(48) 



(49) 



It is worth to mentioning that all relations we obtained in this section restor their respective expressions in the 
previous section when we set Qk — as we expect. Furthermore, one can see that the evolution equation of the GDE 
and the deceleration parameter reduce to the expected formulas in the Einstein's theory [27]. A close look to Fig.© 
reveals some interesting features of the model. First of all, we can see that the EoS parameter, wc in all of the figures 
just can result in an acceleration phase if e < 0.01 which this point is completely consistent with the observational 
evidences [43] . Additionally one can note that in all cases including flat, non flat, interacting and noninteracting 
versions of the GDE model, the EoS parameter of GDE can achieve the phantom cross with suitable choices of e 
and the interaction coupling parameter b. As an example taking = 0.01 favored by observation, fljj = 0.72 for 
the present time, e = 0.002 and b > 0.32 lead to we, < — 1.01. It is worth noting that the existence of the spatial 
curvature in this section does not lead to a considerable difference between the acceleration rate or the deceleration 
rate of the universe with respect to the flat case. However, it should be emphasize that the existence of a non-zero 
curvature in the universe has very important consequences on global topology of the universe and much subtle issues. 



V. DISCUSSION 



Recently, the vacuum energy of the Veneziano ghost field in a time-dependent background was proposed as a kind 
of DE candidate to explain the acceleration of the cosmic expansion. In this model, the energy density of the dark 
energy is proportional to the Hubble parameter H. In this paper we investigated GDE model in the BD framework. 
In this framework we studied the GDE in three distinct cases. At first, we generalized the GDE to BD theory in 
the absence of interaction between DM and GDE in a flat FRW background. In this case the EoS parameter was 
obtained and it was found that for VLd — 0.72, u> — 10000 and e = 0.002 the universe enters a phase of accelerated 
expansion at the late time. Another interesting result from this part is that one can generate a phantom-like EoS from 
a noninteracting GDE in a flat background, which is in contrast to the GDE in standard cosmology. For example, 
taking £1d = 0.72, uj = 10000 and e = 0.007 leads wd — —1.04. In the next two sections we presented the interacting 
GDE in both flat and non-flat BD cosmology. In both of these sections we found that the universe can enter a phase 
of acceleration as well as it can cross the phantom line with suitable choice of constants. For example for = 0.01 
favored by observation, iljj = 0.72 for the present time, e = 0.002 and b > 0.35 lead to Wd < —1.05. For all cases 
we plotted wd versus e and found that crossing from the phantom line is possible provided we take e < 0.01 and in 
the region with e > 0.01 we sometimes even cannot have a phase of acceleration. The important point is that such 
a constraint exist from the observational point of view [43| and we obtained such a condition theoretically. Also we 
investigated the deceleration parameter in different cases. For flat universe in the absence of interaction between DE 
and DM we found that with Qn = 0.72 for present time and choosing e = 0.002 and lj — 10000, we reach q = —0.36 
which shows a consistent value in comparison with observation. Also this parameter have a consistent behavior in all 
cases of interest. 
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